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Abstract. We present the first results of the treatment of grain growth in our 3D, two-fluid 
(gas+dust) SPH code describing protoplanetary disks. We implement a scheme able to reproduce 
the variation of grain sizes caused by a variety of physical processes and test it with the analytical 
expression of grain growth given by Stepinski & Valageas (1997) in simulations of a typical T 
Tauri disk around a one solar mass star. The results are in agreement with a turbulent growing 
process and validate the method. We are now able to simulate the grain growth process in a 
protoplanetary disk given by a more realistic physical description, currently under development. 
We discuss the implications of the combined effect of grain growth and dust vertical settling 
and radial migration on subsequent planetesimal formation. 
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1. Introduction 

Collisions and aggregation govern the first steps of planet formation from micron-sized 
particles to decimetric pre-planetesimals. Observations of protoplanetary disks support 
this mechanism by showing evidence of dust grain growth (e.g. Apai et al. 2004). We 
describe here the implementation of grain growth in our simulations of protoplanetary 
disks. Depending on their relative velocities and material properties (Dominik & Tielens 
1997), collisions between solid particles can make them stick and grow, or conversely 
break them into smaller pieces. 

2. Grain growth in protoplanetary disks 

A first description of the grain growth process is given by Stepinski & Valageas (1997). 
They model a turbulent, vertically isothermal protoplanetary disk, in which gas and dust 
are represented by two separate phases interacting via aerodynamic drag in the Epstein 
regime. Their solid particles are supposed to stick perfectly during collisions, and can 
therefore only grow. The variation of their size s is given by the following analytical 
expression of the growth rate: 
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where Ro is the Rossby number for turbulent motions, a the Shakura & Sunyaev (1973) 
viscosity parameter, p s the density of matter concentrated into solid particles, pd the 
intrinsic density of the grains, C s the local gas sound speed, and Sc the Schmidt number 
of the flow which estimates the effect of gas turbulence on the grains. Sc is defined by 



Sc= (l + fi k t 8 )yi + _3, (2.2) 

where fik is the local keplerian velocity, t s the dust stopping time, v the mean relative 
velocity between gas and dust, and Vt a turbulent velocity. The growth rate depends on 
s through the stopping time 

ts = ^r, (2.3) 

where p g is the gas density. 

The motion of a dust grain and the grain growth process, which both depend on the 
grain size s, are coupled phenomena. The understanding of the global evolution of dust 
in disks therefore requires a numerical treatment. 



3. Grain growth with an SPH code 

Our 3D, bi-fluid, Smooth Particle Hydrodynamics (SPH) code has been developed to 
model vertically isothermal, non self-gravitating protoplanetary disks. Gas and dust are 
treated as two separate phases and are coupled by aerodynamic drag. The code and its 
results on dust migration and settling are presented in Barriere-Fouchet et al. (2005). It 
has also been applied to grain stratification in GG Tau's circumbinary ring (Pinte et al. 
2007) and to gaps opened by planets in the dust phase of protoplanetary disks (Maddison 
et al. 2007; Fouchet et al. 2007; see also talk by Gonzalez, this volume). 

The assumptions we have made are very similar to those of Stepinski & Valageas (1997) 
mentioned in Sect. [21 their prescription for grain growth is therefore easy to implement in 
our code. We allow the grain size s assigned to each SPH particle, assuming it to represent 
the typical size of dust grains at its position, to vary with time following Eq. (|2.ip . We 
take the initial grain size distribution to be uniform. 



4. Results 

We model grain growth in a typical T Tauri disk of mass A/disk = 0.01 M Q , with a 
total dust mass Md us t = 0.01 Mdi s k, around a central star of mass = 1 M . We ran 
simulations with 200,000 SPH particles, evolved over 50,000 years, for a series of initial 
dust grain sizes ranging from sq = 1 /im to 1 mm. Fig. [T] and Fig. [3] show the resulting 
grain sizes in face-on and edge-on views of the disk, for so = 1 /urn. 

We find that grain growth occurs very quickly, especially in the inner disk where the 
density is the highest, as one would expect from Eq. (|2.1[) . In this region, the dust grains 
reach centimetric size (see Fig. [2]) in only a few timesteps. In the outer parts of the disk, 
where the density is far lower, the grains grow much more slowly and their size stays 
below the millimeter. 

An inwards migrating grain coming into a denser region grows almost instantaneously 
to a size characteristic of the grains at its new position, which in turn evolves slowly 
over time. Hence, the global shape of the plot shown in Fig. [2] stays the same as time 
goes on, but the size distribution slowly progresses to larger sizes. The initial size sq has 
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very little influence on the final size distribution, and in particular on the maximum size 
reached. 

Our results are consistent with those obtained by Dullcmond & Dominik (2005). With 
a model solving the coagulation equation in presence of turbulence, they found a very fast 
grain growth in T Tauri protoplanetary disks, depleting very small sizes and producing 
centimetric grains. 

5. Conclusion 

We have implemented in our 3D SPH code a mechanism able to treat grain growth 
and validated it through the use of the simple model of Stepinski & Valageas (1997). 

In accordance with physical intuition, dust grains grow much more quickly in the 
denser, central regions of the disk, where centrimetric sized are reached. The growth 
time is smaller than the migration time and a quasi-stationnary distribution of grain 
size appears in the disk. The small grains are depleted too rapidly to be consistent 
with observations of protoplanetary disks, showing the need to take into account other 
processes such as microscopic interactions between the grains, kinetic energy dissipation, 
porosity, and re-fragmentation. This will be addressed in future work. 
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